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We define a new class of Z' models with neutral flavor-changing interactions at tree level in the 
down-quark sector. They are related in an exact way to elements of the quark mixing matrix due 
to an underlying flavored U(l)^ gauge symmetry, rendering these models particularly predictive. 

The same symmetry implies lepton-flavor non-universal couplings, fully determined by the gauge 
structure of the model. Our models allow to address presently observed deviations from the SM and 
specific correlations among the new physics contributions to the Wilson coefficients Cg jg can be 
tested in 6 —>■ st^l~ transitions. We furthermore predict lepton-universality violations in Z' decays, 
testable at the LHC. 


I. INTRODUCTION 

The Standard Model (SM) of electroweak interac¬ 
tions [1], based on the gauge group Gsm = SU(2)^ x 
U(l)y, has been tested to very high accuracy over the 
last decades. This has strong implications for new 
physics (NP) scenarios, which are then required to in¬ 
volve a very high mass scale, a highly non-trivial fla¬ 
vor structure, or both. Recently the LHCb collaboration 
has performed two measurements involving semileptonic 
b -T si'^t~ transitions which show deviations with re¬ 
spect to the SM expectations: The ratio Rk = Br(i? —>■ 
Kpr)-T Ke^e~) has been measured with a 
central value indicating a violation of lepton universal¬ 
ity at the 25% level [5], implying a 2.6cr deviation, and 
the observable obtained from the angular analysis of 
B — >■ K*pL^pL~ decays differs from the SM expectation 
with 2.9 ct significance [5]. Taking these measurements 
at face value, they require the NP to have the following 
features: 

(i) Sizable contributions to 6 —>■ st^£~ transitions, 

(ii) lepton non-universal couplings. 

Extensions of the SM gauge group by an additional 
U(l)^ factor are among the possible NP scenarios that 
could explain these deviations [IHTU] (other NP interpre¬ 
tations have been discussed e.g. in Refs. [Tlfll4] ). Such 
U(l)^ models have been popular extensions of the SM for 
many years, see Ref. [TS] and references therein. How¬ 
ever, the majority of the extensions considered in the 
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literature are family universal in the quark (and charged- 
lepton) sector, and therefore do not meet condition (ii). 
Models with departures from family universality can give 
rise to large flavor changing Z' interactions, which are 
strongly constrained by current flavor data lillllll]- 
The conditions (i) and (ii) therefore limit considerably 
the viable U(l)^ gauge symmetries. 

Starting a new construction with the minimal particle 
content, i.e. the SM one without right-handed neutrinos, 
the only possible extension is gauging one of the combi¬ 
nations of family-specific lepton number. La — Lp (with 
a,l3 = e,p,T and a ^ (3) [TH]. Since the resulting Z' 
boson couples only to leptons, these models do not meet 
requirement (i) and have been discussed mostly in the 
context of neutrino phenomenology, usually considering 
the — Lr symmetry m- Recent models HIS] cir¬ 
cumvent this problem by introducing additional fermions 
(vector-like quarks). This yields the required couplings 
at an effective level. 

In order to avoid such exotic vector-like quarks, the 
U(l)^ symmetry has to involve both quarks and leptons, 
and should be family-dependent in order to satisfy condi¬ 
tions (i) and (ii). However, such a non-trivial flavor sym¬ 
metry in the quark sector requires the extension of the 
scalar sector in order to accommodate the quark masses 
and mixing angles [501 • been shown in Ref. |5| that 

this option can be realized by including an additional 
complex scalar doublet. This model generates flavor¬ 
changing phenomena in the down-quark sector which 
are related in a first approximation to the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix |^. Flavor¬ 
changing interactions are also present in the up-quark 
sector of this model, but these are not related to CKM 
matrix elements or quark masses. 

We would like our model to have all of its fermion 
couplings related exactly to elements of the CKM ma¬ 
trix. To that aim, we note that there is a class of two- 
Higgs-doublet models (2HDM) that achieves this, known 
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as Branco-Grimus-Lavoura (BGL) models [22]. These 
models have flavor-diagonal interactions in the up-quark 
and charged-lepton sectors, together with flavor-changing 
neutral currents (FCNCs) in the down-quark sector. The 
latter are suppressed by quark masses and/or off-diagonal 
GKM elements in an exact way, thereby providing an al¬ 
ternative solution to the flavor problem in 2HDMs which 
differs radically from the hypothesis of natural flavor con¬ 
servation ( 22 ). 

While in the original BGL model the flavor symme¬ 
try is global, in this work we promote it to a local one. 
This is achieved by charging also the leptons under the 
symmetry, thereby enabling anomaly cancellation. In 
this gauged BGL framework (U( 1 ) 3 ql)i the properties 
of the BGL models are transferred to the gauge boson 
sector: we obtain FGNGs mediated at tree-level by the 
neutral scalar and massive gauge vector bosons of the 
theory, all of which are suppressed by off-diagonal GKM 
elements and/or fermion masses, and therefore naturally 
suppressed. This class of models necessarily exhibits de¬ 
viations from lepton universality due to its gauge struc¬ 
ture. Note that this form of flavor suppression in the 
down-quark sector also appears in certain 3-3-1 mod¬ 
els |7]. However, in our framework this is not a choice, 
but a consequence of the gauge symmetry. 

This paper is organized as foll ows: We formulate the 
U(1 )bgl models in Sec.|^ In Sec. Ill we discuss the con¬ 
straints from flavor and collider data, testing the strong 
correlations present in these models, specifically with re¬ 
spect to the deviations observed by LHCb. We conclude 
in Sec. |IV| The appendices include technical details on 
anomaly cancellation and the scalar sector. 


II. GAUGED BGL SYMMETRY 

A characteristic feature of BGL models |22Q is the 
presence of FGNGs at tree level, which are however 
sufficiently suppressed by combinations of GKM ma¬ 
trix elements. This is a consequence of specific pat¬ 
terns of Yukawa couplings, generated by correspond¬ 
ing charge assignments under a horizontal, family-non- 
universal (BGL-)symmetry. 

The quark Yukawa sector of the model reads {i = 1, 2) 

r, + gi A, 5,4 + h.c., (1) 

where F^ and denote the Yukawa coupling matri¬ 
ces for the down- and up-quark sectors, respectively, 
and = Z(T 2 ‘h*, with the Pauli matrix cr 2 - The neu¬ 
tral components of the Higgs doublets acquire vacuum 
expectation values (vevs) |($°)| = Ui/-\/2. Defining 
V = {vf + measurements of the muon lifetime fix 


V = — 246 GeV. We further define tan/3 = 

V 2 /V 1 and use the following abbreviations: cos/3 = cp, 
sin/3 = s^, tan/3 = tp. 

Choosing an abelian symmetry under which a field ip 
transforms as 

( 2 ) 


the most general symmetry transformations in the quark 
sector yielding the required textures are of the form 

— 2 Xir) + XdR 1 ] , 

4 = diag(A„«,A„fl,Atfl) , (3) 

4 = XdRt, 


with X^R 7 ^ XtR. The Higgs doublets transform as 


=diag(A4>,,A<j,J 

= 2 ^^^® {XuR. — XdR, XtR — XdR) ■ 


( 4 ) 


There are several possible implementations of this sym¬ 
metry, related by permutations in flavor space and ex¬ 
changing up- and down-quark sectors; here, for definite¬ 
ness, the top quark has been singled out, yielding the 
patterns 


* * I 

0 0 oy 


^0 0 0^ 
F2 = 1 0 0 0 


/* * o\ 
Ai = I * * 0 j 
\o 0 oy 


/o 0 o\ 
A2 = 0 0 0 

\o 0 *y 


( 5 ) 


These textures give rise to FCNCs in the down-quark 
sector, which are however suppressed by quark masses 
and off-diagonal elements of the third GKM row [22) . 
This is related to the fact that the GKM mixing ma¬ 
trix V = Ul^iJJdL has a direct correlation in this model 
with the elements of the third row of the rotation of 
the left-handed down quarks, i.e. {UdL) 3 i = V^i. The 
choice of the top quark in Eq. (|^ implies a particularly 
strong suppression of flavor-changing phenomena in light- 
quark systems. Present constraints from AF = 1 and 
AF = 2 quark flavor transitions are accommodated even 
when the scalars of the theory are light, with masses of 
C>(102) GeV |2S|. 

While this model can be implemented through an 
abelian discrete symmetry, this always yields an acci¬ 
dental continuous abelian symmetry which introduces an 
undesired Goldstone boson into the theory [ 22 |- Several 
solutions to this problem have already been proposed in 
the literature ( 22 ] [ 26 ] ; in the following we provide a new 
solution, promoting the BGL symmetry to a local one. 

When gauging a symmetry, special attention must 
be paid to whether it remains anomaly- freej3 BGL 


^ Their extension to include neutrino masses has been discussed in 

Ref. |24| . 2 ^ork we consider only symmetries as “good” if all gauge 
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models are automatically free of QCD anomalies, i.e. 
U(1)^[SU(3)(;;]^ [26]. However, we also need to fulfill the 
anomaly conditions for the following combinations: 

U(1)'[SU(2)J2, U(l/[U(l)y]2, [U(1)TU(1)^, 

[U(l)f, U(l)'[Gravity]2. (6) 

We find that there is no solution for this system within 
the quark sector alone with the charge assignments in 
Eq. Satisfying these anomaly conditions is highly 

non-trivial and requires, in general, additional fermions. 
However, the implementation of the BGL symmetry as 
a local symmetry is possible by adding only the SM lep- 
tonic sector when allowing for lepton-flavor non-universal 
couplings. As in the SM, the cancellation of the gauge 
anomalies then occurs due to a cancellation between 
quark and lepton contributions; the cancellation will 
however involve all three fermion generations, contrary 
to the SM gauge group, for which the anomaly cancella¬ 
tion occurs separately within each fermion generation. 

Gauging the BGL symmetry therefore not only re¬ 
moves the unwanted Goldstone boson, which is “eaten” 
by the Z', but also provides a consistent Z' model with¬ 
out extending the particle content beyond the additional 
scalars, i.e. no vector-like quarks are necessary. The 
flavor structure is extended to the Z' couplings, yield¬ 
ing again strongly suppressed FCNCs in the down-type 
quark sector, while the up-type-quark couplings remain 
diagonal. Especially, since the flavored gauge symmetry 
fixes also the charged-lepton Yukawa matrices to be diag¬ 
onal, we can obtain large deviations from lepton univer¬ 
sality without introducing lepton-flavor violation. Our 
models provide thereby explicit examples for which the 
general arguments given in Refs. muMiiin] do not hold. 

A. Lepton sector 

As emphasized above, the anomaly constraints can be 
fulfilled by assuming that leptons are also charged under 
the U(l)^ symmetry. The resulting charge constraints are 
given in Appendix [A| 

For = 0 fli® mixing between the neutral mas¬ 

sive gauge bosons is suppressed for large tan/3. With 
this choice, we obtain only one possible solution to the 
anomaly conditions up to lepton-flavor permutations, im¬ 
plying the following charge assignments: 

= '^ie = diag(^-^,-^,l^ , 

=diag(^-^,-^,l^ , X[=dmg(^^,^,-3^ , 

anomalies are canceled. It would be also possible to con¬ 
sider anomalous U(l) extensions in the low-energy theory whose 
anomaly graphs get canceled via the so-called Green-Schwarz 
mechanism m- 


= diag Q, y,-3^ , = diag 0^ . 

( 7 ) 

We have normalized all U(l)^ charges by setting XdR = 
1 , thereby fixing also the normalization for the gauge 
coupling g' . In addition to the top, this particular model 
singles out the tau lepton as the only one coupling to $ 2 - 
Permutations of the charges in X[ and A^ give rise to 5 
more models. 

With the charge transformations in Eq. Q, the 
charged-lepton Yukawa sector takes the form 

_^.ie^pto„.^^n,$,eUh.c., (8) 

with 

A 0 o\ /o 0 o\ 

Hi = 0 * 0 , Ha = 0 0 0 , (9) 

\o 0 oy \o 0 

while the quark Yukawa sector remains unchanged. 

We call to attention that no reference to the neutrino 
sector has been made in the previous arguments. In this 
work we are mostly interested in the properties of the 
quark and charged-lepton sectors in the presence of a 
new gauge boson, while we leave a detailed discussion of 
the neutrino sector to future work. For a possible imple¬ 
mentation with the SM fermion content we can build the 
d = 5 Weinberg effective operator which af¬ 

ter spontaneous symmetry breaking induces a Majorana 
mass term for the neutrinos. If this effective operator 
is invariant under the new flavored gauge symmetry in 
Eq. 0, we are led to a — L^. symmetric d = 5 ef¬ 
fective operator and, consequently, Majorana mass term. 
This structure has been studied in Refs. [T9|, and while it 
is not capable of fully accommodating the neutrino data, 
it can serve as a good starting point. 

If we want to improve on this scenario while keeping 
our solution to the anomaly conditions, we may want a 
mechanism that at the effective level already breaks the 
accidental — Lt symmetry. For instance, we may ex¬ 
tend the model with 3 right-handed neutrinos, where one 
of them is not charged while the other two have opposite 
charges. In this way the anomaly solutions remain un¬ 
changed. Choosing the charges of the right-handed neu¬ 
trinos appropriately, the Dirac mass term for the neutri¬ 
nos can be made diagonal and the Majorana mass for the 
right-handed ones — L^. symmetric. Along the lines 
of Ref. |5], coupling an additional scalar to the right- 
handed neutrinos could break the accidental symmetry 
at the effective level. 

We can also envisage other mechanisms that do change 
the anomaly conditions and, therefore, introduce new so¬ 
lutions. 

B. Gauge boson sector 

In order to avoid experimental constraints on a Z' bo¬ 
son, the U(l)^ symmetry must be broken at a relatively 
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high scale, rendering the new gauge boson significantly 
heavier than the SM ones [15]. This can be achieved 
through the introduction of a complex scalar singlet S, 
charged under the BGL symmetry with charge Xs, which 
acquires a vev |(S')| = vsjV^ S> v. The charge of the 
scalar singlet is fixed once the scalar potential is speci¬ 
fied, in our case Xs = —9/8, see Appendix [B] for details. 

As a consequence of the gauge symmetry breaking, 
Gsm X U(1)^ —>■ U(l)gjjj, the neutral massive ^uge bosons 
mix, giving rise to a Lagrangian of the forirj^ 

c=- 

— 1 A — 7^7 — 7 ^ 7 ' 

>7^X1^ 

where and are the SM currents of the photon 
and the Z, respectively. The fermionic piece of the new 
current takes the form 




(xt).PL + (xt).PR 

\ / 27 \ / 27 




3 ’ 


( 11 ) 


with Pl,r = (1 T 75)/2 denoting the usual chiral pro¬ 
jectors, g' the U(l)^ gauge coupling, and the phase 
transformation matrices in the fermion mass basis, i.e. 


'yV _ -y 

/in — /i, 


R ’ 



Xl = Xl, Xl=Xi, and 


9 \vu? v,sv:, v,,v:\ 

+ T VuV:, AtfcV/* 

^ \vtdv:, VtsV.i \vtb\y 


( 12 ) 


The neutral gauge boson mass and mass-mixing pa¬ 
rameters in this model are given by 


tan 2^ = 


2 A 2 


M| - M|, 


(14) 


Since we assume v <C vs, we obtain M|, <C M|/ and 
the resulting mixing angle is small |15j : 

^ A 3. 01“ -1- A o. 01“ / n'nr,n, \ ^ 


e + X,j,^v 2 


9e f g'cpv 


2cwsw 




dtCwSw \ Xlz' 


(15) 

Expanding in this small parameter, we obtain the fol¬ 
lowing leading expression for the hatted weak angle in 
terms of physical quantities 

s 2 „2 _ „2 ,-2 I 

Sw — ~ 72 : 2 “S 


Cw 

In the physical basis the neutral gauge boson masses 
read 




-1 


(16) 


^z,z' — 7 


{( 


M|, + Ml T \ [Ml, - Ml) + 4A4 


and their fermionic currents are given by 

vA* r,(')’/’ p I (')’/’ p 1 / 
ir + iRAiPR Vo, 




^ ~ ^R,ij 

with the couplings (up to O (/)) 


(17) 

(18) 


(A -. 


CwSw 


1 + V 


2 [mi 


- 1 


[Tr-s: 




/ M) 

^X,ij — 


CwSW 


g CwSw cjo/i 
—“ 


C (Ps^ ~ SwQip^ 


(19) 


Ml = 


4s?, 7^2 


Ml,=g'^{vfXl+vlXl) , 
eg' 


A^ = -- 


2cwsw 




(13) 


where e = /A/a is the electromagnetic charge and its ra¬ 
tio with the SU(2)^ coupling g defines 9w = arcsin [e/g), 
which differs from the experimentally measured weak an¬ 
gle due to mixing effects [30]: swcwMz = SwcwMz, 
where Mz is the physical Z mass. We rotate to the 
physical eigenbasis by performing the following orthogo¬ 
nal transformation: 



with O = ( and 

V-'S? 


Here is the third component of weak isospin for the 
left-handed fields (Tjf ” = 0) and the electric charge. 
We also define the vector and axial-vector combinations 
of the Z^''> couplings, 

^V,ij — ^L,ij ^R,ij ’ — ^R/j ^L,ij ’ 

The coefficients ■■ encode corrections to the SM 

Z couplings due to mixing with the Z' which are pro¬ 
portional to g/ at leading order in the mixing an¬ 
gle. Such corrections are restricted to be small {g/ < 
io-'‘) ini ED- In addition to our initial assumption. 
Vs ^ V, which already gives a small / we will assume 
that the mixing angle receives an additional suppression, 
because tan f3 is largej^ This also guarantees that flavor¬ 
changing effects mediated by the Z are negligible com¬ 
pared with those of the Z'. 


® We neglect here a possible mixing coming from the gauge ki¬ 
netic Lagrangian as the mass mixing shows already all relevant 
qualitative effects m- 


^ This limit is quite natural in our model since $2 is coupled to the 
top quark, while $1 is coupled to the light first two generations. 
Therefore, the hierarchy U2 > v\ accommodates well the quark 
mass spectrum. 
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Another important theoretical constraint is obtained 
when requiring the absence of a low-energy Landau pole, 
i.e. an energy scale Alp at which perturbativity is lost. 
This scale can be found from the renormalization-group 
running of the U(l)^ coupling, we have: 

^ = &a'2+0(a'3), (21) 

din 


where the one-loop beta function b contains the charge 
information of the model. It reads 



with / including all fermion degrees of freedom and i = 
1, 2. The Landau pole can be extracted from the pole of 
a'{q^) = g'^(q^)/47r, we obtain 


Alp ~ Mz' exp 


1 

2b a'{Ml,) 


(23) 


For the charges in our model we have b ~ 12.90, see 
Eq. 0. Taking — 4 TeV, we get a Landau pole at 
the Planck scale, App ^ 10^® GeV, for g' < 0.12. Notice 
that this bound is stronger than the naive perturbativity 
bound a'{Ml,) < C>(l/max{|Ari|}®), which in our model 
gives g' < 0.47. We can relax this condition by assum¬ 
ing that additional NP will appear at the see-saw or the 
Grand Unification scale. This way we obtain g' < 0.14 
for Alp > lO^^ GeV and g' < 0.13 for Alp > 10^® GeV. 
The extension of the model to these scales is beyond the 
scope of this work. 


C. Scalar sector 

Here we discuss the main features of the Higgs sec¬ 
tor of the model, a detailed derivation of these results is 
given in Appendix It consists of two complex Higgs 
doublets and a complex scalar gauge singlet. Their prop¬ 
erties are largely determined by the fact that v/vs <C 1. 
The spectrum contains four would-be Goldstone bosons 
giving mass to the electroweak (EW) gauge vector bosons 
{Mw,Mz,Mz'}, a charged Higgs and four neutral 
scalars (three CP-even {Hi,H 2 ,H^} and a CP-odd A). 

We obtain a decoupling scenario with a light SM- 
like Higgs boson, ~ 0{v'^), three heavy quasi¬ 

degenerate states coming mainly from the scalar dou¬ 
blets ~ Ml ~ ^H± C>(u|) and another heavy 

CP-even Higgs coming mainly from the scalar singlet 
Ml, - o(^^|)g The Yukawa interactions of the scalar 


^ It is also possible to motivate scenarios that avoid decoupling 
based on an enhanced Poincare symmetry protecting the weak 
scale m- These can give rise to a rich scalar sector at the EW 
scale and a pseudo-Goldstone boson in the spectrum. 


bosons reflect this decoupling. They can be written as 


-CyD pi dR + huyr m + y/ en] 

ip^HkA 

M _ 

H- - uJiNldL + (-LNien) + h.c. 

(24) 

Neglecting corrections of 0{v’^/vD we obtain 


- RL 

^ f ’ 

J y 

yA _ 

“ V ’ 


\aH2 _ _ " ' J 

^ f 

J y 


Nf 


— 7 




(25) 


with (/ = u, d, £). Here fermions have been rotated to the 
mass eigenbasis. The diagonal matrices Df^u,d,i contain 
the fermion masses, and the matrices Nf are given as 


= ? (Dd),^ - f? + y) ivA{v)s,{Dd)j,, 

•' Vi \Vi V2 J 

Nu = — diag(TO„, rric, 0) - —diag(0, 0, rrit ), 

Vi V2 

Vo V^ 

Ne = —diag(TOe, m^, 0)-diag(0,0, nir) ■ (26) 

Vi V2 


The couplings of the scalar to the fermions are 
additionally suppressed, since they are generated only 
through mixing with the doublet degrees of freedom: 

~ 0{v/vs)- The matrix is non-diagonal in flavor 
space, giving again rise to tree-level FCNCs in the down- 
quark sector controlled by the CKM matrix V. Further¬ 
more, as shown in Appendix]^ the constraints from the 
U(l)^ symmetry render the scalar potential CP invariant. 
CP violation in our model is therefore uniquely controlled 
by the CKM phase, like in the SMj^ 


D. Model variations 

As mentioned in section [II A[ we have six model imple¬ 
mentations according to the different lepton flavor per¬ 
mutations of the U(l)^ charges, with identical quark and 
scalar charges given in Eq. 0. These permutations are 
performed in the basis where the lepton mass matrix is 
diagonal and the eigenvalues are correctly ordered. This 
way each permutation corresponds to the simultaneous 
interchange of the charges of both left- and right-handed 
leptons. It is straightforward to check that the permu¬ 
tation of just one set of charges, while satisfying the 
anomaly conditions as well, would give rise to a non¬ 
diagonal mass matrix and reduce to one of the six cases 
considered here after diagonalization. 


® We do not consider strong CP violation in this work. 
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In order to present the predictions for each of these 
models, we introduce the generic lepton charges 


XiL = 

9 

4’ 

X 2 L = 

21 

T ’ 

XsL = 

XlR = 

9 

2 ’ 

X 2 R = 

15 

y ’ 

X 3 R = 


such that each implementation is labeled by (e, r) = 
{i, j, k). For instance, the model presented in Eq. 0 is 
now labeled as (1, 2, 3). 

The additional possibilities to implement the quark 
and neutrino sectors have been discussed in Secs. HI] 
an d jirai respectively. 

III. DISCUSSION 

In this section we discuss the phenomenology of the 
gauged BGL models introduced in the previous section. 
The phenomenology of a scalar sector with a flavor struc¬ 
ture identical to the one of our model has been analyzed 
in Ref. Since we are assuming a decoupling scalar 

sector, we can naturally accommodate a SM-like Higgs 
at 125 GeV; the heavy scalars are not expected to yield 
sizable contributions to flavor observables in general. We 
therefore focus on the phenomenological implications of 
the Z' boson for this class of models. In the cases where 
the constraints depend on the choice of lepton charges 
in our model, we give here the most conservative bound 
and discuss the differences in Secs. IIII Cl and HITTO 

A. Low-energy constraints 

Despite the large mass of several TeV, the Z' bo¬ 
son yields potentially significant contributions on fla¬ 
vor observables, due to its flavor-violating couplings 
in the down-type quark sector. However, because of 
M^/Mf, <0.1% and the GKM suppression of the flavor¬ 
changing Z' couplings, these contributions can only be 
relevant when the corresponding SM amplitude has a 
strong suppression in addition to Gp- This is specifically 
the case for meson mixing amplitudes and electroweak 
penguin processes which we will discuss below. Further 
examples are differences of observables that are small in 
the SM, for example due to lepton universality or isospin 
symmetry. 

The particular flavor structure of the model implies 
a strong hierarchy for the size of different flavor transi¬ 
tions: \V*Vtd\ ~ A5 « \V;^Vtb\ ~ A3 « \V,lVtt\ ~ A2 
(where A ~ 0.226). However, since the SM amplitudes 
often show a similar hierarchy, the relative size of the Z' 
contribution has to be determined on a case-by-case ba¬ 
sis. We obtain for example similar bounds from the mass 
differences Arrid^s in the Bd,s systems and e^. Here we 
include exemplarily the constraint from Bg mixing, while 
leaving a detailed phenomenological analysis for future 
work. 


Given the high experimental precision of Arug [34] , the 
strength of the corresponding constraint depends com¬ 
pletely on our capability to predict the SM value. The 
limiting factors here are our knowledge of the hadronic 
quantity '/Bs (for which substantial progress is ex¬ 
pected soon, see Ref. [3S]) and the relevant GKM ele¬ 
ments. Using the corresponding values from Refs. |36j 
and [37], respectively, we conclude that contributions to 
Anis up to 20% remain possible at 95% confidence level 
(GL); from this we obtain the limit Mz'jg' ^ 16 TeV. 
Tree-level scalar contributions have been neglected, since 
their effect cancels to a very good approximation in the 
decoupling limit considered here |38) . 

Regarding the input for the GKM parameters, we make 
the following observations in our models: (i) Charged- 
current tree-level processes receive negligible contribu¬ 
tions. (ii) The mixing phase in B^^a mixing and the ratio 
Ams/Arud remain SM-like. (iii) The additional direct 
CP violation in H —>■ TnT,pTT,pp and B J/'hK is neg¬ 
ligible. These observations imply that in our context the 
global fits from Refs. [371133] remain valid to good ap¬ 
proximation. 

In our models the Z' boson couples to muons and 
will contribute to neutrino trident production (NTP), 
—>■ vNp.'^pT. We expect dominance of the Z' vec¬ 
tor current for our models. The NTP cross section nor¬ 
malized to the SM one was calculated in Refs. HiaEn]. 
Using a combination of the latest NTP cross-section mea¬ 
surements we obtain a bound on the parameter combina¬ 
tion Mz'lg' [H], which is however weaker than the one 
from Ba mixing. 

Atomic parity violation (APV) measurements also 
place bounds on additional neutral gauge bosons coupling 
to electrons and light quarks. In particular, the precise 
measurement of the weak charge from APV exper¬ 

iments is in a reasonable agreement with the prediction 
of the SM and can be used to set bounds on the Z' bo¬ 
son of our models j33]. The bound we obtain on Mz' jg' 
from the latest determination of the -v^gak charge 

is however again weaker than the limit from Ba mixing. 

Further potentially strong bounds stem from electric 
dipole moments (EDMs) and the anomalous magnetic 
moment of the muon. EDMs from the scalar sector oc¬ 
cur at the two-loop level; while this alone does not render 
them necessarily sufficiently small, the additional sup¬ 
pression and cancellations in the decoupling limit imply 
very small effects |3^|33]. Anomalous magnetic moments 
can be present at the one loop level, but are again too 
small in 2HDMs near the decoupling limit jSS] El] • Gon- 
cerning the new gauge boson, one-loop contributions to 
EDMs require different phases for left- and right-handed 
Z' couplings to quarks and/or leptons [33]; this feature 
is not present in our model, therefore the Z' contribu¬ 
tions are again at the two-loop level and sufficiently sup¬ 
pressed, given the large Z' mass. In the case of anoma¬ 
lous magnetic moments we have non-vanishing one-loop 
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contributions, which for the muon read |46j 


a 


NP 


9' 


47r2 Ml, 



(28) 


where and = {X^r - 

X^r)/2. Using the bound on Mz'/g' from mix¬ 

ing, we get < 1.3 x 10“^^ for all charge assignments, 
which is smaller than the current theory uncertainty in 
the prediction of [3T] . 


B. Direct searches 


The obvious way to search directly for a Z' is via a 
resonance peak in the invariant-mass distribution of its 
decay products. At the LHC this experimental analysis 
is usually performed by the ATLAS [47] and CMS [48] 
collaborations for Z' production in the s-channel in a 
rather model-independent way, but assuming validity of 
the narrow-width approximation (NWA), negligible con¬ 
tributions of interference with the SM [49], and flavor- 
universal Z' couplings to quarks. Under these assump¬ 
tions, the cross section for pp —)• Z'X —)• ffX takes the 
simplihed form isniEii 


48s L 


clwu (s, M|,) -f c{wd (s, M|,) 


(29) 


where the functions Wu,d are hadronic structure factors 
that encode the information of the Drell-Yan production 
processes of the Z' and their QCD corrections (for a pre¬ 
cise definition of these functions we refer the reader to 
Ref. [ST]). The model-dependent part of the cross sec¬ 
tion is contained in the coefficients Cu,d'- 


cl^g'^{Xl^ + Xl^)BT{Z'^ ff) , 
ci^9'\Xl^+Xl^)BT{Z'^ff) . 


While the assumptions for these expressions are not ex¬ 
actly fulfilled in our model, they are applicable when ne¬ 
glecting the small contributions proportional to the off- 
diagonal CKM matrix elements. The reason is that under 
this approximation the Z' couplings to the first two gen¬ 
erations - which yield the dominant contribution to the 
Drell-Yan production of the Z' - are universal and fla¬ 
vor conserving in the quark sector, see Eqs. Q and (12). 
Note that X^l — XdL up to small corrections propor¬ 
tional to the off-diagonal CKM matrix elements that we 
are neglecting}^ 

For high values of g' the ratio Tz>/Mz' can be quite 
large, spoiling the NWA. The constraint is therefore only 


^ Flavor changing Z' couplings between the first two quark gen¬ 
erations are suppressed by lUdLj* | ~ A® while violations of uni¬ 
versality are suppressed by lUdl^ ~ lUaP ~ A'*, with 

A ~ 0.226. 


applicable for g' < 0.2, such that Tz' jMz’ does not ex¬ 
ceed 15%. However, for the mass range accessible so far 
such a large coupling is typically excluded by the con¬ 
straint Rk^ as will be shown below. Combined exclusion 
limits on the Z' at 95% CL using the dimuon and dielec¬ 
tron channels have been provided by the CMS collabo¬ 
ration in the (c„,Cd) plane [IHj. Using these bounds to¬ 
gether with the constraint from Rk, we find that a Z' bo¬ 
son is excluded in our models for Mz' < 3-4 TeV, depend¬ 
ing on the lepton charge assignments, while g' should be 
0(10“^). These bounds are stronger than those from 
many Z' benchmark models, for which current LHC data 
typically give a lower bound of Mz' > 2-3 TeV (27] HHj . 
The reason is a combination of a sizable value for g' jMz' 
required by the constraint from Rr and the combination 
of Z' couplings to fermions appearing in our model, see 
Eq. 

When the Z' is too heavy to be produced at a given 
collider, M\, ^ s, indirect searches can be carried 
out by looking for effects from contact interactions, 
i.e. effective operators generated by integrating out 
the heavy Z'. At the LHC, the corresponding expres¬ 
sions hold for Mz' ^ 4.5 TeV, and limits can be ex¬ 
tracted e.g. from searches for the contact interactions 
{q"i^iPL,Rg){il^ Pl,r() [52|. These are presented in terms 
of benchmark scenarios considering a single operator with 
a specific chirality structure. From the 95% CL limit pro¬ 
vided in Ref. [52] for the operator {qjij^PLq){e'y^P rc) we 
extract Mz'jg' ^ 22 TeV. Searches for contact interac¬ 
tions of the type {e'-^fiPL,Re){Jl^PL,Rf) were performed 
at LEP [5S]. These are sensitive to a Z'-boson coupling 
to electrons assuming Mz' ^ 210 GeV and the result¬ 
ing limits are again given for benchmark scenarios. In 
this case we extract Mz' jg' > 16 TeV from the 95% CL 
bound on the operator (e 7 “PRe)(/ 27 aPRp,). It is impor¬ 
tant to note that the extracted bounds from contact in¬ 
teractions using benchmark scenarios do not capture the 
full dynamics of our models once the Z' boson is inte¬ 
grated out. Interference effects due to operators with 
different chirality structures will be relevant in general, 
a detailed analysis of these effects is however beyond the 
scope of this work. These bounds should therefore be 
taken as a rough guide to the sensitivity of contact inter¬ 
actions searches at LEP and the LHC to the Z' of our 
models. However, in general they do not put relevant 
bounds as long as g' is not too large. 


C. Discriminating the different models 

The constructed class of models has two important fea¬ 
tures which imply specific deviations from the SM: con¬ 
trolled FCNCs in the down-quark sector and violations 
of lepton universality. Thanks to the very specific fla¬ 
vor structure, the two are strongly correlated in each of 
our models. While the flavor-changing couplings are uni¬ 
versal in the models we discuss here, the lepton charges 
differ. Correspondingly the models can be discriminated 








by processes involving leptons, specifically rare leptonic 
and semileptonic decays of B mesons, which can test both 
features. 

The effective Hamiltonian describing these b —>■ st^t~ 
transitions reads 


- ^ ^ {c/ol + 0o'f) + 

^ i 


(31) 

where 


O 9 = (ire), 

= {s^i^PRb) (eri ), 

Olo = irnPLb) (irid) , 

0(0 = irnPRb) (jrid) , 

= mtisPRb){ee ), 

o'l = mi,{sPLb){ee ), 

Op = TOb(sPp6)(/75€) , 

0 '^ = mi,{sPLb){hd) ■ 


(32) 


The piece is the effective Hamiltonian for the ra¬ 

diative transitions including dipole operators m- New 
physics contributions to arise in our models from 

one-loop diagrams mediated by Z' and Higgs bosons. 
Both of them are very small due to the assumed mass 
scale of several TeV m- 

At the 6-quark scale, the SM contribution to the re¬ 
maining part reads C™ ~ “C'lo^ — 4.2 V£, while all 
other contributions are negligible. The chirality-flipped 
operators Og receive negligible contributions also in 
our models since = 0, while all others receive con¬ 
tributions from Z' or Higgs exchanges [I2ll38]- The Z' 
contribution to Cg j^g is given by 


'NP^ ^ 

TTv"^ 

J d J e 
^L,sb^V,U 

9 — 

aV^gVti, 

Ml 

'NP^ ^ 

■nv^ 

J d Me 

10 — 

c^VfgVti, 

Ml 


(33) 


where Cf = In Table|^we show the correla¬ 

tions between the Wilson coefficients C^^fg in our models 

and provide as a function of g' jMz' by introducing 

a conveniently normalized parameter Kg: 


cr- s < X 10< (1^)' = 4 X 605 TeV» (^-0 . 

(34) 

This allows the direct estimation of C^fg in terms of 
g'/M z 4 useful for phenomenological purposes. 

The CP-even Higgs H 2 and the CP-odd Higgs A can 
give sizable contributions to Og and Op, respectively. 
For muons we hav^ 


O5 ~ ~^P 


-2TT{t/3 + 

aMjj 


(35) 


The correlation Cs = —Cp was expected given the assumed 
decoupling in the scalar sector m- 


Table I. Correlations among the NP contributions to the ef¬ 
fective operators Og^io- 


Model 

^10 /^9 

/^NPe /y^lNP/x 
O 9 /O 9 

/^NPe //-yJNP/x 
^10 7^9 

^9 

(1,2,3) 

3/17 

9/17 

3/17 

-1.235 

(1,3,2) 

0 

-9/8 

-3/8 

0.581 

(2,1,3) 

1/3 

17/9 

1/3 

-0.654 

(2,3,1) 

0 

-17/8 

-3/8 

0.581 

(3,1,2) 

1/3 

-8/9 

0 

-0.654 

(3,2,1) 

3/17 

-8/17 

0 

-1.235 


Here we have denoted the quasi-degenerate masses of H 2 
and A by Mp = ~ M^, see Appendix]^ for details 

of the scalar sector. The coupling is again model- 

dependent: for models (1,2,3), (2,1,3), 

(3,1, 2) and (3, 2,1), while = -0mp for (1, 3, 2) 

and (2, 3,1). The suppression by the muon mass renders 
these contributions negligible, apart from observables in 
which also the SM and Z' contributions receive this sup¬ 
pression, e.g. Bs —)■ Higgs contributions to O'0 

will be additionally suppressed by a factor mg/mi, com¬ 
pared to the corresponding non-primed operators and are 
neglected in the following. 

The angular distributions of neutral-current semilep¬ 
tonic 6 —?► s£^i~ transitions allow for testing these coef¬ 
ficients in global fits. Furthermore, they provide precise 
tests of lepton universality when considering ratios of the 
type 


Rm = 


Br{B -» Mg+g, ) sM 
Br(H^Me+e-) 


1 + 0{ml/ml), 


(36) 


with M G {K, K*, Xg, Ko(1430),...} [5S], where many 
sources of uncertainties cancel when integrating over 
identical phase-space regions. 

Additional sensitivity to the dynamics of NP can be 
obtained via double-ratios [55] . 


(37) 

It was shown in Ref. [SS] that the dependence on the 
NP coupling to left-handed quarks cancels out in Rm- 
Since Cg^ig = 0 in our models we have Rk* = Rx„ = 

77^0 ( 1430 ) = 1) providing an important test of the flavor 
structure of our models. 

The possible hadronic final states for these ratios yield 
sensitivity to different NP structures, thereby providing 
complementary information. The discriminating power 
of these observables has been shown recently in Ref. [T^ 
within the framework of leptoquark models and in more 
general contexts in Refs. [laiMiisz]. 

To analyze the deviations from flavor-universality we 
define Rm — 1 -I- A -1- S, where E is the pure contribution 
from NP and A the one from the interference with the 
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SM. These quantities are given by 


A = 






Re 


{^SM 


+ Re{cS“(c'foP^)*}-(Ai^e) 


(38) 


S = 


a 


NP/i 


a 


NP^t 

10 


lei'll 


Ir'SMi 

r^io 




and are valid to a very good approximation given the 
present experimental uncertainties. 

Our models also give clean predictions for the rare de¬ 
cay modes B —>■ {K, K*, Xs}^^ |17) . Tree-level Z'-boson 
contributions to these decays are generically expected of 
the same size as in 6 —> s££ transitions. However, the 
modes with neutrino final states do not distinguish be¬ 
tween our different models, since there is no sensitivity to 
the neutrino species. We obtain again a universal value 
for all ratios = Bt{B —>■ Mvv)/Br{B —>■ MPi^)|sm, 
due to the fact that = 0. For the same reason the 
average of the K* longitudinal polarization fraction in 
B —)■ K*vv is not affected by the Z' exchange contribu¬ 
tion m- The enhancement of turns out to be rel¬ 
atively small, 0(10%) for g' ~ 0.1 and Mz' ^ O(TeV). 
The reason is again the sum over the different neutrino 
species: due to the different charges, the Z' contribution 
interferes both constructively and destructively, leaving a 
net effect that is smaller than in the modes with charged 
leptons. 

The leptonic decays > £~^£~ constitute another 

sensitive probe of small NP effects. Within the SM, these 
decays arise again at the loop level and are helicity sup¬ 
pressed. Due to the leptonic final state, the theoretical 
prediction of these processes is very clean. They receive 
Z'- and Higgs-mediated contributions at tree-level in our 
model through the operators Ofg g p [13] , 


Br{Bs ) 

Br(Hs ^ /r+At-)™ 


1 - 0.24C'iT'^ - 


where ~ 7.7mb- For the range of model parameters 
relevant here, the Z' contribution to Bs —> is very 

small with respect to the SM, ~ 1% for Mz' ~ 5 TeV 
and g' ^ 0.1. Larger contributions are possible due to 
scalar mediation for the models (1, 2, 3), (2,1,3), (3,1, 2) 
and (3,2,1), given that Cs — —Cp ex in these 
cases. Taking Mp ~ 10 TeV and tp ^ 30 for exam¬ 
ple, one obtains a suppression of Br(i3s —)■ gXgT) of 
about 10% relative to the SM. In our model the ratio 
Bi[Bbi —>■ g~)/BT{Bs ) remains unchanged 

with respect to the SM to a very good approximation. 

If a Z' boson is discovered during the next runs of the 
LHC [58|, its decays to leptons can be used to discrimi¬ 
nate the models presented here. We dehne the ratios 


^ Z' ^ //) 

^Z'^ f'f) ’ 


( 40 ) 


Table II. Model-dependent predictions for the ratio gi/p. 


Model 


Me/r 


(1,2,3) 

45/149 

45/32 

149/32 

(1,3,2) 

45/32 

45/149 

32/149 

(2,1.3) 

149/45 

149/32 

45/32 

(2,3,1) 

149/32 

149/45 

32/45 

(3.1,2) 

32/45 

32/149 

45/149 

(3,2,1) 

32/149 

32/45 

149/45 


where again the dominant sources of uncertainty can¬ 
cel, rendering them particularly useful to test lepton 
universality. This possibility has also been discussed in 
Refs. [331 . In our models we have 


fJ-b/t — 






Xfip + Xpp 


(41) 


The hrst ratio is fixed in our models, gb/t — 1, while 
large deviations from gi/p = 1 are possible. The model- 
specific predictions for the ratio gi/p are shown in Ta¬ 
ble |n| We note that there is a correlation between Rp 
and ge/fi'- an enhancement of Rk (as present in models 
(1,3,2), (2,1,3) and (2,3,1)) implies an enhancement 
of ge/fj. and the other way around. Violations of lepton 
universality can therefore be tested by complementary 
measurements of flavor observables at LHCb or Belle H 
and Z' properties to be measured at ATLAS and CMS. 


D. Interpretation of 6 s£^£ anomalies 

The LHCb collaboration recently performed two mea¬ 
surements of decays with b —>■ s£'^£~ transitions that 
show a tension with respect to the SM expectations. 
The first measurement is that of the ratio Rk introduced 
in the last section, where LHCb measures m ijLHCb ^ 
0 .745to;o74 ± 0.036 for G [1,6] GeV^, a 2.6 cr deviation 
with respect to the SM prediction [55l|60l|6T|. The second 
measurement is the angular analysis of H —>■ K* g'^g~ 
decays [3|, where an excess of 2.9cr is observed for the 
angular observable in the bin G [4,6] GeV^ which 
has been constructed to reduce the influence of form fac¬ 
tor uncertainties [5^ . 

Attributing the measurement of Rk solely to NP, i.e. 
assuming that it will be confirmed with higher signifi¬ 
cance in the next run of the LHC and/or by additional 
Rm measurements, it is the first way to exclude some of 
our models, since it requires sizable non-universal contri¬ 
butions with a specific sign. The flavor structure in each 


® The bin £ [6, 8] GeV^ shows nominally a tension with identi¬ 
cal significance; however, this bin is considered less theoretically 
clean, due to the larger influence of charm resonances. 
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of our models is fixed, so half of them cannot accom¬ 
modate Rk < 1, namely (1, 3, 2), (2,1,3) and (2, 3, 1). 
This is illustrated in Fig.[2 where it is additionally seen 
that a large deviation from Rk = 1 , as indicated by the 
present central value and Icr interval, can actually only 
be explained in two of the remaining models. This strong 
impact shows the importance of further measurements of 
Rm ratios. 


Table III. Model-dependent bound on from the Rk mea¬ 

surement. The constraint from Bs mixing is taken into ac¬ 
count. 


Model 


Cr^{2a) 

(1,2,3) 

(3,1,2) 

(3,2,1) 

[-0.93, -0.43] 
[-1.20,-0.53] 

[-2.92,-0.61] 

[-1.16,-0.17] 

[-1.54,-0.20] 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

SfV/WzlTeV-'] 


Figure 1. Model-dependent predictions for Rk as a func¬ 
tion of g' /Mz'. The recent measurement of Rk by the LHCb 
collaboration is shown at Icr and 2a. Constraints from Bs 
mixing are also shown at 95% CL. 


In Fig. we show the constraints from the Rk mea¬ 
surement for the remaining models (1,2,3), (3,1,2) and 
(3, 2,1). The allowed regions are consistent with the con¬ 
straint from Sg-meson mixing. LHC searches for a Z' 
boson exclude values of Mz> below 3-4 TeV, as discussed 
in Sec. IIIB the corresponding areas are shown in gray. 
We also show the theoretical perturbativity bounds ob¬ 
tained from the requirement that the Landau pole for 
the U(l)^ gauge coupling appears beyond the see-saw or 
the Grand Unification scales, i.e. Alp > 10^^ GeV and 
Alp > 10^® GeV, respectively. 

Regarding the angular analysis in 5 —>■ K* fr^ pT, the 
situation is more complicated. The wealth of informa¬ 
tion provided in this and related measurements requires 
a global analysis, which is beyond the scope of this work. 
However, a number of model-independent studies have 
been carried out [81 I63U67] . These analyses differ in terms 
of the statistical methods used, treatment of hadronic 
uncertainties (see also Ref. [SH]), and consequently in the 
significance they find for the NP hypothesis over the SM 


.NP/i 


-1 


one. However, they agree that a contribution Cg 
can fit the data for without conflicting with other ob¬ 
servables. Additional contributions from e.g. can 


be present, but are less significant and tend to be smaller. 
When using Table [I| to translate the Rk meas urem ent 


into a bound on Cg ^ in our models, see Table 


HI 


the 


ranges are perfectly compatible with the values obtained 
from P 5 , as also observed for other Z' models [U EU E?] . 
This is highly non-trivial given the strong correlations in 
our models and will allow for decisive tests with addi¬ 
tional data. 

Let us further comment on the implications of the re¬ 
cent measurements of H))decays [69], 


Br(H, 


Br(Hs 
Br(i?d ^ 

Br{Bd /r+^-)SM 


./-lexp 

T' ) _ n '7C-I-0.20 

+ ' 0 - 0.18 > 


(42) 


= 3.7 


- 1 - 1.6 
- 1 . 4 ' 


Both results seem to hint at a deviation from the SM, 
however in opposite directions. A confirmation of this sit¬ 
uation with higher significance would rule out our models 
which predict these ratios to be equal; however, the un¬ 
certainties in the mode are still large. Regarding the 
Bs mode, which is measured consistent with the SM pre¬ 
diction, it should be noted that the result depends sen¬ 
sitively on the value adopted for |I 4 h|j where the value 
from inclusive B -T Xciv decays was chosen in the SM 
calculation |7D|. In any case, as discussed in Sec. 
a potential shift could be explained in the context of our 
models by scalar contributions, which are, however, not 
directly related to the Z’ contributions discussed above. 

The models that accommodate the b —>• s£^£~ data 
can be further discriminated using the observables de¬ 
scribed in the last section, notably using the ratios 
provided a Z' boson is discovered at the LHC. 


me 


IV. CONCLUSIONS 

The class of family-non-universal Z' models presented 
in this article exhibits FCNCs at tree level that are in 
accordance with available flavor constraints while still 
inducing potentially sizable effects in various processes, 
testable at existing and future colliders. This is achieved 
by gauging the specific (BGL-)symmetry structure, in¬ 
troduced in Ref. [22] for the first time, which renders the 
resulting models highly predictive. 

The particle content of the models is minimal in the 
sense that the extension of the U(l)^ symmetry to the lep- 
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0.20r 


0.15 


0.00 


0.20 


0.15 


g'0.1 


0.00 


2 3 4 5 6 7 

Mz'fTeV] 


ALP>10^^GeV 
ALP>10’®GeV 


Model (123) 

2 (T 


CMS (95% CL) - 
7TeV - 
8TeV - 



0.20 


0.15 


g'0.10 


Model (321) 
□ 2-r 

■ 

CMS (95% CL) 
7TeV 
■ 8TeV 


2 3 4 5 6 7 

Mz'[TeV] 


Figure 2. Regions allowed at la and 2a by the Rk measurement in the {Mz'^g'} plane for the models (1,2,3), (3,1,2) and 
(3,2,1). Exclusion limits from Z' searches at the LHC are shown in gray. The black lines indicate bounds from perturbativity 

of g'- 


ton sector allows to restrict the fermion content to the 
SM one. The only additional particles are then a second 
Higgs doublet, a scalar singlet and the Z' boson, all of 
which are heavy after spontaneous symmetry breaking, 
due to the large mass scale for the singlet. The anomaly 
conditions largely determine the charge assignments un¬ 
der the U(l)^ symmetry; the remaining freedom is used 
for two phenomenologically motivated choices, leaving 
only six possible models which are related by permuta¬ 
tions in the lepton sector. 

The main phenomenological features of these models 
can be summarized as follows: 


Any significant deviation from the SM in an observable 
with Z' contributions allows for predicting all other ob¬ 
servables, for instance the values for obtained from 

Rk, see Tables |T] and [ml Further characteristic predic¬ 
tions in our models include the following: 


• The absence of flavor-changing Z' couplings to 
right-handed quarks implies Rm = !> i-e. all ra¬ 
tios defined in Eq. (36) are expected to be equal, 
Rk = Rk* = Rx^ = ... The same holds for the 
ratios va B ^ Mvv decays. 


1. FCNCs at tree level in the down-quark sector, 
which are mediated by heavy Z' gauge bosons and 
neutral scalars, controlled by combinations of CKM 
matrix elements and/or fermion mass factors. 

2. Non-universal lepton couplings determined by the 
charges under the additional U(l)^ symmetry. 

3. No FCNCs in the charged-lepton or up-quark sec¬ 
tors. 

4. Complete determination of the U(l)^ sector up to 
two real parameters, Mz' and g', where all observ¬ 
ables at the electroweak scale and below depend 
only on the combination g' jMz'- 

Present data already strongly restrict the possible pa¬ 
rameter ranges in our models: direct searches exclude 
Z' masses below 3 — 4 TeV and the constraint from B 
mixing implies Mz' jg' > 16 TeV (95% CL). Theoretical 
bounds from perturbativity give an upper bound on the 
value of the gauge coupling, e.g. g' < 0.14 for a Landau 
pole beyond the see-saw scale. Nevertheless, three of our 
models can explain the deviations from SM expectations 
in 6 —>■ t~ transitions seen in LHCb measurements 

mia, while the other three are excluded (at 95% CL) by 
Rk < 1. These findings are illustrated in Figs. and 


• If a Z' is discovered at the LHC, measurements of 
(7{pp Z' ^ £ili)/a{pp Z' ^ fj^j) can be 
used to discriminate between our models, due to 
the specific patterns of lepton-non-universality, see 
Table ni 

• Leptonic down-quark FCNC decays are sensitive 

to the Higgs couplings in our model, specifically 
Bd,s —t ■ Double ratios of Bg and B^ decays 

are again expected to equal unity. 

In the near future we will therefore be able to dif¬ 
ferentiate our new class of models from other Z' mod¬ 
els as well as its different realizations from each other. 
This will be possible due to a combination of direct 
searches/measurements at the LHC and high-precision 
measurements at low energies, e.g. from Belle H and 
LHCb. Further progress can come directly from theory, 
e.g. by more precise predictions for ^ or . 

As a final remark, we recall that in this minimal imple¬ 
mentation of gauged BGL symmetry neutrinos are mass¬ 
less. Additional mechanisms to explain neutrinos masses 
can introduce new solutions to the anomaly equations; 
these new variants are subject to future work. 
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Appendix A: Anomaly cancellation conditions 

In this appendix we present the restrictions on the 
charges derived from the requisite gauge anomaly can¬ 
cellation for quarks transforming according to Eq. @ 
and allowing the leptons to transform in the most gen¬ 
eral way: 


= diag (VeL, X^L,XrL) , 

Xf, = diag{XeR,X^R,XrR) . 

As mentioned in Sec.[^ the QCD anomaly condition is 
automatically satisfied within BGL models. On the other 
hand, the anomaly conditions involving a single U(l)^ 
gauge boson together with SU(2)j^ or U(l)^ ones, and the 
mixed gravitational-U(l)^ anomaly are non-trivial and 
read 


3 9 

-^ 221 ' = SXuR + ^XtR + -XdR + , 




Alll' — — -XuR — -^XtR — -^XdR 




otR 


(A2) 


AgGV = ^ {‘2‘XaL — X, 


aR) 


ct—e,^,T 


The triangle diagrams involving two or three U(l)^ gauge 
bosons result in more complicated conditions for the BGL 
charges, which include a quadratic and a cubic equation: 


Aii'i' = 


- 2^1r - ^^tR + -^^dR + XuRXdR 


-^XtRXdR — ^ 


{xL - xln) , 


Av 


-\[^xIr- 


XtR + iX, 


dR 


-XdR { 2 X 1 ^ + Xj^) - Xl^ { 2 X^r + X^r)] 

+ ^ {2X1^ - Xl^) . (A3) 

Satisfying these conditions with rational charges is non¬ 
trivial and there is only one class of solutions giving an 
anomaly-free model, up to lepton flavor permutations: 

XuR = —XdR — -Xf^R , XtR = —AXdR + -X^R , 

1 5 

XeL = XdR + ^Xfj^R , Xfj_L = —XdR + -^X^R , 

9 

X^L — '^XdR X^R , 

1 4 

XeR = 2XdR + -X^R , XrR = IXdR — -X^R . 

(A4) 

At this point we have two free charges in the above rela¬ 
tions, i.e. XdR and X^^r. Using = 0, as discussed 
in Sec. |TTj yields one relation between the two charges. 
The remaining free charge, e.g. XdR, is fixed by some 
normalization convention. 


Appendix B: Scalar potential 


The Higgs doublets are parametrized as 



( ^ 

+ Pi + "iVi)J 


(* = 1 , 2 ). 

(HI) 


Their neutral components acquire vevs given by ($°) = 
e’‘^iVjjy/2 with Vi > 0 and (vf + = v. In full 

generality we can rotate away the phase of 4>i, leaving 
the second doublet with the phase Q = 82 — Oi. We 
parametrize the complex scalar singlet as 


S = 


-ias/2 (**S + Ro+ Up) 

V 2 


(B2) 


with Vs > 0. 

Since we have in our model three scalar fields, 
{$i,$2 ,5'}, the phase-blind part of the scalar potential 
has a U(l)^ global invariance. In order to avoid massless 
Goldstone bosons, we need to charge the S field in such 
a way that the phase-sensitive part breaks this symme¬ 
try down to U(l)^ = U(I)y X U(l)^ Gauge-invariant 
combinations can be built only from ${$2 and S, leav¬ 
ing us with the possibilities $|$ 2 S', 4)|$25'^ and com¬ 
binations involving complex conjugates. Eor concrete¬ 
ness we choose to be invariant, which imposes 

Xs = —9/8. The scalar potential then reads 


U=m?|$,|2 + ^|$,|4 + A3|$inci>2|2 

+ A4|4>I$2p + (B3) 

- ^ ($I$252 + $^2 4>i(5*)2) . 
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All the parameters of the scalar potential but ^3 are real 
due to hermiticity. The parameter ^3 has been chosen 
to be real and positive by rephasing the scalar held S 
appropriately. 

The vacuum expectation value of the potential is given 
by 


r 


+ 2Ai4 + ^4 + (*^3 + A4)c|^ 

+ 2A25^ + 2 v^ + f^lC^ + 623 ^^ + d 2 V^ 

— 2 63 Cpsi 3 V^ cos{as — 9) . (B4) 


1/v <C 1. Including the leading corrections in 1/v one 
obtains 


7^s 


with 


/ cp 

-S0 

\ ^23^13 — 


S/3 Wi3 \ 

C/3 OJ 23 , 

— (t^23C/3 + W13S/3) 1 J 

(B9) 


Wi3 = — 


2 dic^ — 63 S 213 + 2(52 s| 
2vd2 


^354/3 + 2(di — ^2)52/3 
2v{2d2S2/3 - S3) 


(BIO) 


Here we have dehned the ratio v = vg/v. The stability ^he resulting masses for the CP-even scalars are given 
of the vacuum requires that 


By convention we choose ag = 9, the other possibili¬ 
ties in ag = 9 mod 27r simply amount to an unphysical 
rephasing of the scalar held S. Stability also requires 
dVo/dvi = 0, to which the only non-trivial solution reads 


^2 


ml 


— [—diC^ -|- (53C/3S/3 — 623 '^ — d2'5^] ) 

[—4 (A3 -l- A4) s^j — 4 Aic^ -|- 53tgv^ — 2 i 5 i'!}^] , 

— —4 (A3-l-A4) — 4A2 s|- l-— 2(52'0^ . 

(B6) 




4s2/3 


Mh - ^ 


(Bll) 


with 


A — {2d2Xi — 51)0^ + {Sic^ + (52S^)(53S2/3 

2i5id2 + <5i 2 (r2 \ 3 4 

*2/3 ~ (^2 ~ 292 X 2 ) 31 : 


(B12) 


+ d2{X3 -|- A4)s2^ . 


The exact expression for their Yukawa couplings in 
Eq. pll) is 


vYj: ’’ = [ci3{TZg)kl + 3i3{TZg)k2] Df 

+ [3/3{'R-g)kl — Ci3{TZg)k2] Nf ■ 


(B13) 


The scalar potential is then determined in terms of 10 
unknown parameters {vg, /3, Ai_4, d 2 , di_3}, since 9 = ag 
does not appear explicitly. 

The masses of the physical CP-odd boson A and the 
charged scalar are given by 


Ml = 


SsVg 

4 



3^\ 

{)2 J 


SsVg 

4 s 2/3 




4(1 / ^ _ 2 A 4 
4 V*2/3 v"^ 


Ssvg 
4s2^ ’ 


(B7) 


respectively. The physical states Hi,2,3 are expressed in 
terms of {pi, p 2 , Ro} via an orthogonal transformation, 


Finally, the vacuum solution in our models allow for 
complex vevs but that is not sufficient to have sponta¬ 
neous CP violation. In the weak gauge sector CP viola¬ 
tion is manifest through the invariant Ti[Hd, Hy]^ [71j . 
In the CP-invariant scenario, i.e. with real Yukawa ma¬ 
trices, the Hermitian combinations Hy^d are given by 

2Hu = yfAiAj + VIA 2 A 2 + uiU2(AiA^ -|- A2Af )ce 
-k iuii;2(AiA^ - A2Al)se , 

2Hd = u^PiPf + VIT 2 TI + viV2{T2Tl + rir^)ce 

-k iViV2{T2^'l - PiP^ss • 

(B14) 



(B8) 


The mass matrix for the CP-even bosons can be diago¬ 
nalized analytically in a perturbative expansion around 


In our model we get Tr [H^/, H„]^ = 0, implying the ab¬ 
sence of CP violation in the gauge interactions when the 
only phase is carried by the scalar vev. As a consequence, 
the source of CP violation for the weak currents in our 
model is present in the Yukawa couplings and will appear 
in the observables through the CKM mechanism. 
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